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1 Introduction 

This note focuses on some issues for the analysis of a system of degenerate cross-diffusion 
partial differential equations (PDEs). This family of models are encountered in a wide 
variety of contexts, such as population dynamics [5], biology mm, chemistry or materials 
science. The application we have in mind here is the modeling of the evolution of the 
concentration of chemical species composing a crystalline solid. The functions, that are 
the solutions of the system of PDEs of interest, represent the local densities of the different 
components of the material, and thus should be nonnegative, bounded and satisfy some 
volumic constraints which will be made precise later in the note. These systems are 
useful for instance for the prediction of the chemical composition of thin solid films 
grown by Chemical Vapor Deposition (CVD) [ 8 ]. In this process, a solid wafer is exposed 
to gaseous precursors, corresponding to the different species entering the composition 
of the film, which react or decompose on the substrate surface to produce the desired 
deposit. This process generally occurs at high temperature and takes several hours, so 
that the diffusion of the different atomic species within the bulk of the solid has to be 
taken into account in addition to the evolution of the surface of the film. 

More precisely, let d E N* (the space dimension), T > 0 (the duration of the pro¬ 
cess), and for all t E [0,r], let C denote a bounded open domain occupied by 
the solid at time t. Let us assume that there are n -|- 1 different atomic species (with 
n E N*) composing the material, whose local concentrations are respectively denoted by 
Uo{t, x),ui{t, x), • • • , Unit, x) for t E (0, T) and x E D*. 

The system of equations we consider here reads: 






where for all 0 < z 7 ^ j < n, the non-negative real number Kij = Kji > 0 denotes 
the cross-diffusion coefficient of atoms of type i with atoms of type j. This system of 
equations can be formally derived as the continuous limit of some discrete stochastic 
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lattice hopping model. We do not explain the derivation of these equations here for the 
sake of brievity and refer the reader to [DEI SI for more details. 

The initial condition {uq,Ui, ■ ■ ■ ,u^) E L^(no; of this system is assumed from 

now on to satisfy the constraints 

n 

y 0 < i < n, for almost all x E fioj ^ 0; (^) = 1 Ui{0,x) = u^{x). 

i=0 

( 2 ) 

The last constraint = 1 is a natural volumic constraint which encodes the 

fact that each site of the crystalline lattice of the solid has to be occupied (vacancies 
being treated as a particular type of atomic species). 

From a physical point of view, it is naturally expected that the solutions (uq, • • • , Un) 
of such a system satisfy similar constraints as the initial condition, i.e. for all t E (0, T), 

n 

V 0 < z < n, Vf E [0,T], for almost all x E Ui{t,x) > 0 and E Ui{t,x) = l. (3) 

i=0 

Indeed, the last equality can be infered at least formally from ([T|) since by summing all 
these equations for 0 < z < n, we obtain that dt '“*) ~ 0 - 

Despite their importance in chemistry or biology, it appears that the mathematical 
analysis of these systems, taking into account constraints ([3]), is quite recent miiniiiDE]. 
Let us first consider the case of a fixed domain, i.e. when D* = Dq for rth t E [0,T], and 
of no flux boundary conditions 

Kij{ujVui — uiS/uj) • n = 0 on (0, T) X V 0 < z < n, (4) 

where n denotes the outward normal unit vector to When all the coefficients Kij 

are identically equal to some constant iL > 0 , it can be easily seen that m boils down 
to a set of decoupled heat equations 

V 0 < z < n, dfUi — diVa; [KVxUi] = 0 in (0, T) x Dq, 

with Neumann boundary conditions KS/xUi • n = 0 on (0,T) x 5Do whose analysis 
becomes trivial. However, at least to our knowledge, the first proof of existence of 
global weak solutions of © satisfying constraints m with non-identical cross-diffusion 
coefficients is given in | 2 ] for n = 2 with coefficients K^j such that KiQ > 0 for z = 1 , 2 
and Ki 2 = K 21 = 0. These results were later extended in [B] to a general number of 
species n E N* with cross-diffusion coefficients satisfying K^q > 0 and Kij = 0 for all 
1 E f 7 ^ J < ri; the authors proved in addition the uniqueness of such weak solutions. 
In |5], the case n = 2 with arbitrary positive coefficients > 0 is covered, though no 
uniqueness result is provided. The main difficulty of the mathematical analysis of such 
equations relies in the bounds Q , which are not obvious since no maximum principle can 
be proved for these systems in general. In all the articles mentioned above, the analysis 
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framework used by the authors is the so-called boundedness by entropy method, which was 
first used in |21 in the particular case mentioned above and then further developped by 
Jiingel in [I] to a more general theoretical setting. The main idea is to write the above 
system of equations as a formal gradient flow and derive estimates on the solutions 
(mQ) ■ ■ ■ using the decay of some well-chosen entropy functional. The diffusivity 
matrix [9j obtained for these systems is not a concave function of the densities though, 
so that standard gradient flow theory arguments (such as the minimizing movement 
method) cannot be applied in this context. In the first part of this note, we prove that 
for any n G N*, in the case when Kij = Kji > 0 for all 0 < i 7 ^ j < n, the above system 
can be analyzed using the general theoretical framework introduced in |3]. The proof 
heavily relies on ideas of [5]. Uniqueness of the solutions remains an open issue. 

Interestingly, it seems that very few works deal with the mathematical analysis of such 
cross-diffusion systems in the case of moving boundary domains, despite their importance 
for instance in the modelisation of the CVD processes mentioned above. In the second 
part of this note, we present preliminary existence results for one-dimensional systems 
with moving boundaries, extending existence results obtained in the general theoretical 
framework introduced by Jiingel [1]. 

2 Existence on a fixed domain: formal gradient flow struc¬ 
ture 

For the sake of completeness, let us recall here Theorem 2 of jl]. 

Theorem 2.1 (Theorem 2 of [1]). Let D be an open domain of ML such that D C (a, b)^ 
for some a, 6 G M. Let A : u £ D ■= (^ij(w))i<i,j<n £ be a matrix-valued 

functional defined on D satisfying A G and the following assumptions: 

(HI) There exists a bounded from below convex function h G C^(I1,M) such that its 
derivative Dh : D —)• is invertible on R"'; 

(H2) For all 1 < i < n, there exist a* > 0 and > 0 such that for all z = 

(zi, • • • , G R"" and u = (tti, • • • ,UnY' £ D, 

n 

D‘^h{u)A{u)z > ai{ui)‘^zf, 

i=l 

where either ai{ui) = a*{ui — 0 )™'*“^ or ai{ui) = af{h — 

(H3) There exists a* > 0 such that for all u £ D and 1 < i, j < n such that nij > 1, it 
holds that \Aij(u)\ < a*\aj{uj)\. 

Let £ L^{LIq‘,D). Then there exists a weak solution u with initial condition to 

dtu = divx{A{u)Vxu) in (0,T) x Oq, and A{u)Vu • n = 0 on (0,T) x OLIq, (5) 

such that for {t,x) £ {0,T) x Oq; u{t,x) £ D with u £ L^((0, T); (Oq, R”)) and 
5inGL2((0,T);FF-Htlo;M-)). 
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In this first section, we prove that the system of cross-diffusion equations introduced 
above, on a fixed domain Qq, 


dtUi - div^ 




0<j^i<n 


= 0 in (0, T) X rio) V 0 < z < n, ( 6 ) 


and with no flux boundary conditions (H]), can be analyzed using the general theoretical 
framework introduced in |3]. Actually, the following existence result is a direct conse¬ 
quence of Theorem 2 of [1]: 


Corollary 2.2. Let us assume that for all 0 < i ^ j < n, Kij = Kji > 0 and 
let := (tig, • • • ,M^) G L^(no; satisfying eondition Then, there exists 

a weak solution u = (zio, • • • ,zin) to with initial condition vf' satisfying u G 

L2((0,r);i?i(f7o;M’^+i)), dtu e and 

n 

V 0 < z < n, for almost all {t,x) G [0, T] x Oq, Ui{t,x) > 0 and E Ui{t,x) = 1. 

i=0 


Indeed, if (zzq, • • • , Un) is a solution to ([ 6 ]) with boundary conditions ([4]), it holds that 
dt ~ ^ almost all {t,x) G (0,T) x Qq. Assuming that the initial condition 

(zzq, • • • , zz°) satisfies m, the above system can then be reformulated equivalently, writing 
zzq := 1 — a function of the concentrations (zzi, • • • ,zz„) only as follows: for 

all 1 < z < n. 


dtUi — div 




+-^*o((l -/o)VzzjzziVp) = 0 in (0,T) x Oq, 
~ UiVuj) + KiQ{{l - p)Vui + UiVp) • n = 0 on (0,T) X Sflo, 

(7) 

where p := Ui = l-UQ. 


In our case, D := {(ui,--- ,Un) G < 1} (0,1)”- For any u = 

(ui, • • • ,Un) G D, we denote by p := 'Ya=i 'T'he system ([7|) introduced above can be 
rewritten under the form ([5]) using the notation u := (zzi, • • • ,Un)'^, with A : u & D 
A{u) := {Aij{u))i<ij<n G defined as follows: for all 1 < z,j < n. 


-^n('a) — ^ ^ i^ij Kio^Uj -|- KiQ and Ajj(zz) — if z ^ j. 


Let us prove that A satisfies the assumptions of Theorem 12.11 with the following 
entropy functional, which is the same as in umai: 

n 

h : u & D h{u) := Ui Inzzj -|- (1 — /?) ln(l — p). ( 8 ) 

i=l 
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The function h G C^{D;'R) nC^(D;]R) (thus is bounded on D), is strictly convex 
on D, and its derivative Dh : D —)• M”' is invertible. For all u £ D, Dh{u) = (Inuj — 
ln(l — p))i<i<n and for all w G M”, Dh~^{w) = ( . ^ ) . As a consequence, 

h satisfies assumption (HI) of Theorem 12.11 

Let us now prove that assumption (H2) of Theorem 12.11 is satisfied with ^ for 

all 1 < f < n. Thus, there will be no need to check assumption (H3) for the existence 
result to hold and Corollary 12.21 will be a direct consequence of Theorem 12.11 To this 
aim, we follow the same strategy of proof as the one used in [5]. We are going to prove 
that if the coefficients Kij are assumed to be strictly positive then, there exists a > 0 
such that for all ti G H, 


H{u)A{u) > aA(tt), where H{u) 


D‘^h{u), A{u) := diag 



and a := min Kij. 
0<i^j<n 


(9) 

This inequality implies (H2) with rrii = ^ and a* = and ai{u) = for all 

1 < f < n. Let u £ D. We have for all 1 < i,j < n, 

Hii{u) = — + and if i ^ j. 

Ui 1- p 1-/9 

Introducing P{u) := the matrix such that for all 1 < i,j < n, 


Pii{u) = l-Ui and Pij{u) = -Ui if i / j. 


it holds that H{u)P{u) = A(n). Thus, we can write H{v)A{u) — aK{u) = H{u){A{u) — 
aP{u)). We can easily check that A{u) — aP{u) = A{u) + aD{u), where A{u) has 
the same structure as A(u) but with diffusion coefficients Kjj — a instead of K^j, and 
D{u) := {Dij{u))i<ij<n where Dij{u) = Uj for all 1 < i < n. 

On the one hand, H{u)D{u) = where Z := {Zij)i<cij<n with Zjj = 1 for all 

^ < i, j ^ n. Since the matrix Z is a positive matrix, so is H{u)D{u). 

On the other hand, since h is strictly convex on D, H{u)A[u) is positive if and only 
if M{u) := A{u)H{u)~^ is positive. Indeed, for all z G M"", we have H{u)A{u)z = 
{H{u)z)'^ (^A{u)H{u)~^^ {H{u)z). Besides, it can be easily checked that M(m) = (M^-(«))!<*j<„, 
where for all 1 < < n, 

Mii{u) = {KiQ-a){l-p)ui+ ^ {Kij-a)uiUj and Mij{u) = -{Kij-a)uiUj if j / i, 

and that this matrix is indeed a positive matrix. Hence we have proved inequality Q, 
assumption (H2) and Corollary 12.21 


3 A simple one-dimensional model with moving boundary 

As announced in the introduction, the second part of this note deals with a model 
describing the evolution of the concentration of different chemical species in a solid thin 
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film during a CVD process. In this case, the evolution of the boundary of the film has 
to be taken into account as well as the diffusion phenomena occuring in the bulk. We 
introduce and analyze a simplified one-dimensional model, considered as a preliminary 
step before tackling more complicated systems. 

In this section, d = 1 and at time t = 0, the domain occupied by the solid is an 
interval Qq = (0, cq) for some initial thickness cq > 0. The external fluxes of volatile 
precursors during the deposition process are denoted by ((/>o, • • • , (j)n) £ L°°((0, T); 
and are assumed to satisfy > 0 for all t G (0,T) and 0 < z < n. At time t G (0,T), 
the domain occupied by the solid is equal to = ( 0 ,e(t)) where the depth e{t) of the 
film is assumed to obey an evolution law characterized only by the values of the external 
fluxes: 

e{t) = eo + / 

■^0 *=o 

The concentrations (uq, • • • , Un) satisfy the evolution law given in ([1]) with the following 
boundary conditions: for all 0 < z < n and t G (0,T), 


I ^ KijiujOxUi UidxUj) I (t)0) 

\0<jj^i<n J 


0 and 


Kij{ujdxUi-UidxUj) {t,e{t))+e'{t)ui{t,e{t)) 


These boundary conditions ensure that the total mass of each chemical species is con¬ 
served in the sense that for all t G (0, T), dt Ui(t, x) dx^ = 

Actually, we can prove the following existence result in a slightly more general frame¬ 
work, directly inspired from the setting of [1]. 

Proposition 3.1. Let D := {(zzi,--- ,zZn)^ G < 1} ^ (0,1)” and 

A : u G D A{u) := (Ajj(zz))i<jj<„ G he a matrix-valued functional defined on 

D satisfying A G and assumptions (H1)-(H2)-(H3) of Theorem \2. 1\ with an 

entropy functional h satisfying in addition 

(H4) h G L^(D). 

Let eo > 0, G L^((0, cq); T>) and G L°°((0, T); We denote by 

(j) := {(fi, • • • , (fn)'^■ Then there exists a weak solution u with initial condition u^ to 

' e{t) = eo-\- fg Er=o 

dtu - dx(A(u)dxu) = 0 for t e (0,T) and X e (0,e(t)), . . 

{A{u)dxu){t, e{t)) + e'{t)u{t, e{t)) = ^{t) for t G (0, T), 

^ {A{u)dxu){t,0)=0fort£{0,T). 


such that for {t,x) G (0,r) x (0,e(t)), u{t,x) G H with u £ L‘^{{0,T); H^{{0,e{t))]W^)) 
anddtu G L2((0,T);F-i((0,e(t));M”)). 


Note that the main difference between the assumptions of Proposition 13.II and Theo- 
rem l2.1l is the fact that the entropy function h is asked in addition to belong to L°°{D] M). 
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Indeed, a crucial point in our proof is that this functional h has to be bounded from above 
and from below on D, which is indeed the case for the function h defined in ([8]). We refer 
the reader to [T] for a proof of this result. The idea of the proof is to define the problem 
onto a reference fixed domain and derive entropy estimates which enable to carry on the 
same analysis as in |3]. 

4 Conclusion 

In this note, we have presented preliminary existence results for a cross-diffusion model 
with a moving boundary domain, modeling the evolution of the concentrations of different 
chemical species in a thin film solid layer grown by CVD. 

In the first section, we have slightly extended results of [5] in the sense that we have 
proved that the system of PDEs ([6]) with boundary conditions (H]) defined on a fixed 
domain, with coefficients Kij = Kji > 0 for all 0 < i ^ j < n, fell into the analysis 
framework developped by Jiingel in [1]. In the second section, we have proposed and 
analyzed a simplified one-dimensional model in order to describe the growth of a thin 
him layer during a CVD process, as well as the diffusion of the different chemical species 
inside the bulk of the solid. The uniqueness of the solution of such a system of equations 
remains an open issue in both cases. 

These results should be seen as preliminary steps before tackling more challenging 
systems of PDEs. In particular, the analysis of more realistic two-dimensional and three- 
dimensional models for CVD deposition taking into account surface diffusion phenomena 
is currently work in progress. 
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